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The absorption spectrum of V838 Mon in 2002 February - 
March. I. Atmospheric parameters and iron abundance/ 
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ABSTRACT 

We present a determination of the effective temperatures, iron abundances, and mi- 
croturbulent velocities for the pseudophotosphcrc of V838 Mon on 2002 February 25, 
and March 2 and 26. Physical parameters of the line forming region were obtained in 
the framework of a self-consistent approach, using fits of synthetic spectra to observed 
spectra in the wavelength range 5500-6700 A. We obtained T eff = 5330 ± 300 K, 5540 
± 270 K and 4960 ± 190 K, for February 25, March 2, and March 26, respectively 
The iron abundance log iV(Fe)=— 4.7 does not appear to change in the atmosphere of 
V838 Mon from February 25 to March 26, 2002. 
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1 INTRODUCTION 

The peculiar variable star V838 Mon was discovered during 
an outburst in the beginning of 2002 January (Brown 2002). 
Two further outbursts were then observed in 2002 February 
(Munari et al. 2002a; Kimeswenger et al. 2002; Crause et 
al.2003) and in general the optical brightness in V-band of 
the star increased by 9 mag. Since 2002 March, a gradual 
fall in V-magnitude began which, by 2003 January, was re- 
duced by 8 mag. The suspected progenitor of V838 Mon was 
identified by Munari et al. (2002a) as a 15 mag F-star on 
the main sequence. Possibly V838 Mon might have a B3V 
companion (Desidera & Munari 2002), but it could be a 
background star. The discovery of a light echo (Henden et 
al. 2002) allowed an estimate of the distance to V838 Mon 
and, according to recent works based on HST data (Bond 
et al. 2003; Tylenda 2004) its distance is 5-6 kpc. If these 
estimations are correct, at the time of maximum brightness 
V838 Mon was the most luminous star in our Galaxy. 

Details of the spectral evolution of the star are described 
in Kolev et al. 2002; Wisniewski et al. 2003; Osiwala et 
al. 2002). During outbursts (except for the last) the spec- 
trum displayed numerous emission lines with P Cyg pro- 
files, formed in the expanding shell and around an F- or 
A-star (Kolev et al. 2002). On the other hand, absorption 
spectra appropriate to a red giant or supergiant were ob- 
served in quiescent periods. Strong lines of hydrogen, D lines 
of sodium, triplets of calcium and other elements show P 
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Cyg profiles. They have similar profiles and velocities vary- 
ing from —500 km s _1 in late January to —280 km s 
in late March (Munari et al. 2002a). Since the middle of 
2002 March, the emissions are considerably weakened and 
the spectrum of V838 Mon evolved to later spectral classes. 
In middle of 2002 April, there were present some lines of 
TiO; in May the spectrum evolved to the "very cold" M- 
giant (Banerjee & Ashok 2002). In October Evans et al. 
(2003) characterized it as a L-supergiant. 

Recently Kipper et al.(2004) found for iron group ele- 
ments [m/H] = —0.4, while abundances of lithium and some 
s-process elements are clearly enhanced. This results was 
obtain using the static LTE model. 

These results are very dependent on the model atmo- 
sphere and spectrum synthesis assumptions. 

The nature of the outbursts remains a mystery. Possible 
explanations include various thermonuclear processes (very 
slow nova, flare post-AGB), and the collision of two stars 
(Soker & Tylenda 2003). Munari et al. (2002a) suggested 
that V838 Mon is a new type of a variable star, because 
comparison with the closely analogous V4334 Sgr and M31 
RV has shown significant enough differences in the observed 
parameters. 

In this paper we discuss the results of the determina- 
tion of iron abundance and atmospheric parameters of V838 
Mon. These we obtained from an analysis of absorption spec- 
tra of V838 Mon on 2002 February 25 and March 2 and 26. 
The complexity and uniqueness of the observed character- 
istics of V838 Mon practically excluded a definition of the 
parameters of the atmosphere using conventional methods, 
based on calibration on photometric indices, ionization bal- 
ance, profiles of hydrogen lines. Indeed, the presence around 
the star of a dust shell, and the uncertain determination of 
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interstellar reddening (from Eb-v =—0.25 to Eb-v =—0.8 
Munari et al. 2002a), affects the U — B and B — V colours. 
Emission in the hydrogen lines provides severe problems for 
their application in the estimation of effective temperature. 
Moreover, both the macroturbulent motions and expansion 
of the pseudophotosphere merges the numerous lines in wide 
blends. As a result, a single unblended line in the spectrum of 
V838 Mon cannot be found at all, and any analysis based on 
measurements of equivalent widths is completely excluded. 

The observational data used in this paper are described 
in section 2. Section 3 explains some background to our work 
and some details of the procedure used. We attempt to de- 
termine T e fl , the microturbulent velocity Vt and the iron 
abundance log JV(Fe) in the atmosphere of V838 Mon in the 
framework of the self-consistent approach in section 4. Some 
results are discussed in section 5. 
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2 OBSERVATIONS 

Spectra of V838 Mon were obtained on 2002 February 25 
and March 26 with the Echelle+CCD spectrograph on the 
1.82m telescope operated by Osservatorio Astronomico di 
Padova on Mount Ekar (Asiago), and freely available to 
the community from http:/ /ulisse. pd.astro.it/V838Mon/ A 
2 arcsec slit was used with fixed E-W orientation, produc- 
ing a PSF with a FWHM of 1.75 pixels, corresponding to 
a resolving power close to 20000. The detector was a UV 
coated Thompson CCD 1024 x 1024 pixel, 19 micron square 
size, covering in one exposure the wavelength range 4500 to 
9480A (echelle orders #49 to #24). The short wavelength 
limit is set by a 2 mm OG 455 long-pass filter, inserted in the 
optical train to cut the second order from the cross-disperser. 
The wavelength range is covered without gaps between ad- 
jacent echelle orders up to 7300 A. The spectra have been 
extracted and calibrated using IRAF software running un- 
der Linux operating system. The spectra are sky-subtracted 
and flat-fielded. The wavelength solution was derived simul- 
taneously for all 26 echelle orders, with an average r.m.s of 
0.18 km s^ 1 . The 8480-8750 A wavelength range of these Asi- 
ago spectra has been described in Munari et al. (2002a,b). 

Another set of spectra (R ~ 32000) for March 2 was 
obtained with the echelle fibre-fed spectrograph on the 1.9- 
m SAAO telescope kindly provided for us by Dr. Lisa Crause 
(see Crause et al. 2003 for details). 



3 PROCEDURE 

To carry out our analysis of V838 Mon we used the spectral 
synthesis techniques. Our synthetic spectra were computed 
in the framework of the classical approach: LTE, plane- 
parallel media, no sinks and sources of energy inside the 
atmosphere, and transfer of energy provided by the radia- 
tion field and by convection. 

Strictly speaking, none of these assumptions is 100% 
valid in atmosphere of V838 Mon. Clearly we have non- 
static atmosphere which may well have shock waves mov- 
ing trough it. Still we assumed that in any moment the 
structure of model atmosphere of V838 Mon is similar to 
model atmospheres of supergiants. Indeed, temporal changes 
of the absorption spectra on the days were rather marginal. 



Figure 1. Comparison of expansion (V^xp = 160 km s ), rota- 
tional (v * sini = 80 km s ) and macroturbulent (V ma cro = 50 
km s -1 ) profiles used in this paper to convolve synthetic spectra. 

Most probably, for this object, we see only a pseudophoto- 
sphere, which is the outermost part of an expanding enve- 
lope. Therefore, our first goal was to determine whether it 
is possible to fit our synthetic spectra to the observed V838 
Mon spectra. 

At the time of the observations the spectral class of 
V838 Mon was determined as K-type (Kolev et al. 2002). 
Absorption lines in spectrum of V838 Mon form compara- 
tively broad blends. 

Generally speaking, there may be a number of broad- 
ening mechanisms: 

• Microturbulence, which is formed by small scale (i.e 
r< 1) motions in the atmosphere. In the case of a super- 
giant, Vt usualy does not exceed 10 km s _1 . In our analysis 
we determined Vt from a comparison of observed and com- 
puted spectra. 

• Stellar rotation. Our analysis shows that, in the case 
of V838, we should adopt v * sini = 80 km s _1 to fit the 
observed spectra. This value is too high for the later stages 
of stellar evolution, for obvious reasons. In reality rotation 
cannot contribute much to the broadening of lines observed 
in spectra of most supergiants. 

• Expansion of the pseudophotosphere of the star. Asym- 
metrical profiles of expansion broadening can be described, 
to a first approximation, by the formula 



G(v, A, AA) = const * AA^l - (AA/A * V C x P /c) 2 , (1) 

where I4x P is the expansion velocity. The observed spectra 
can be fitted with V cxp = 160 km s _1 (see Fig. 0. 

From February 25 to March 2 the brightness of the star 
was approximately constant, but at the time of the obser- 
vation on March 26 there was a small drop in visual bright- 
ness (see Kolev et al. 2002) . In Fig. |2]we show a comparison 
of the observed spectra of V838 Mon at the times of the 
increase in luminosity (February 3; Kipper & Klochkova, 
private communication) and at approximately constant lu- 
minosity (February 25, March 2 and March 26). They differ 
drastically - the spectrum of February 3 is severely veiled by 
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Figure 2. Spectra of V838 Mon observed on February 3, Febru- 
ary 25, March 2 and March 26 2002 
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Figure 3. Comparison of the spectrum of V838 Mon and that 
of Arcturus, convolved with macroturbulent profile Vmacro = 50 
km s 



emission: many lines are observed in emission. This demon- 
strates that effects of the radial expansion of the line-forming 
layers were not significant for the dates of our data and for- 
mally obtained value V oxp = 160 km s _1 is not real. 

• Macroturbulence. After the large increase of luminos- 
ity in 2002 January-February, large scale (i.e. of magnitude 
t > 1) macroturbulent motions should be very common in 
the disturbed atmosphere of V838 Mon. Our numerical ex- 
periments showed that, to get appropriate fits to the ob- 
served spectra taking into account only macroturbulent ve- 
locities, we should adopt Vmacro ~ 50 km s _1 . 

In any case, for the times of our observations the spectra 
of V838 Mon resemble the spectra of "conventional" super- 
giants. Our V838 Mon spectra for February 25, March 2 
and 26 agree, at least qualitatively, with the spectrum of 
Arcturus (K2 III), convolved with macroturbulent velocity 
profile, given by a gaussian of half-width Vmacro = 50 km 
s" 1 (Fig.0. 

The observed emissions in the cores of the strongest 
lines are formed far outside, perhaps at the outer boundary 
of the expanding envelope, i.e. in the region which is heated 
by shock wave dissipation. 

As result of our first numerical experiments, we con- 
cluded that the spectra of V838 Mon in 2002 February - 
March were similar to the spectrum of a normal late (su- 
per)giant, broadened by strong macroturbulence motions 
and/or expansion of its pseudophotosphere. Unfortunately 
we cannot, from the observed spectra, distinguish between 
broadening due to the macroturbulence and expansion (see 
next section). 

It is worth noting that the observed spectra of V838 
Mon are formed in a medium with decreasing temperature to 
the outside, i.e. in the local co-moving system of co-ordinates 
the atmosphere, to a first approximation, can be described 
by a "normal" model, at least in the region of formation of 
weak or intermediate strength atomic lines. 



3.1 Fits to observed spectra 

We computed a sample of LTE synthetic spectra for a grid of 
Kurucz (1993) model atmospheres with T e fl = 4000 - 6000 
K using the WITA612 program (Pavlenko 1997). Synthetic 
spectra were computed with wavelength step 0.02 A, micro- 
turbulent velocities 2-18 km s _1 with a step 1 km s _1 , 
iron abundances logTV(Fe) = —5.6 — * — SSdex 1 , with a step 
0.1 dex. Then, due to the high luminosity of the star, we 
formally adopt log g = 0. Synthetic spectra were computed 
using the VALD (Kupka et al. 1999) line list. For atomic 
lines the line broadening constants were taken from VALD 
or computed following Unsold (1955). 

For the dates of our observations lines of neutral iron 
dominate in the spectra. Fortunately, they show rather weak 
gravity /pressure dependence, therefore the uncertainty in 
the choice of log g will not be important in determining our 
main results; the dependence of the computed spectra on 
T e fi is more significant (see Fig. 0J. The computed syn- 
thetic spectra were convolved with different profiles, and 
then fitted to the observed spectra following the numeri- 
cal scheme described in Jones et al. (2002) and Pavlenko & 
Jones (2002). 

In order to determine the best fit parameters, we com- 
pared the observed residual fluxes r^ bs with computed values 
H{ h +l r - We let Hf s = f F* c ° r *G(y)*dy, where Ff- eof is the 
theoretical flux and G(y) is the broadening profile. In our 
case G(y) may be wavelength dependent. To get the best fit 
we find the minima per point of the 2D function 

S(f e , / g ) = S(l - H synt /H° bs f. 

We calculated these minimization parameters for our 
grid of synthetic spectra to determine a set of parameters / a 
(wavelength shift parameter) and / g (convolution parame- 
ter). 

The theoretical spectra were convolved with a gaussian 
profile. Our convolution profile is formed by both expan- 
sion and macroturbulent motions. We cannot distinguish 
between them in our spectra. To get a numerical estimate 

1 in the paper we use the abundance scale iVj = 1 
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T pff =4000 K 
T B)f =5000 K logg=0 
T pff =5000 K logg=1 
T eff =6000 K 



Table 1. Parameters of the solar atmosphere 
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Figure 4. Dependence of computed spectra on T e fj and log g 



of the broadening processes in the pseudophotosphere, we 
use a formal parameter V g , which describes the cumulative 
effect of broadening/expansion motions. 

The parameters / a and / g were determined by the min- 
imization procedure; the procedure was carried out for dif- 
ferent spectral regions. We selected for analysis 6 spectral 
orders in the interval 5600-6700A. In the red, spectral lines 
are blended by telluric spectra, and are of lower S/N. In 
the blue the blending of the spectra are rather high. Our 
main intention was to obtain a self-consistent solution sep- 
arately for different echelle orders, and then compare them. 
If we could obtain similar parameters from different spec- 
tral regions it can be evidence of the reality of the obtained 
solution. 
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Figure 5. The best fit to the solar spectra. 



4 RESULTS 
4.1 The Sun 

To be confident in our procedure, we carried out a similar 
analysis for the Sun. For this case we know the solar abun- 
dances and other basic parameters, therefore our analysis 
provides an independent estimation of the quality of our 
procedure: 

• From the solar atlas of Kurucz et al (1984) we ex- 
tract spectral regions corresponding to our observed orders 
of V838 Mon; 

• we convolve the solar spectra with a gaussian of 
Vmacro = 50 km s -1 . 

• we carried out a spectral analysis of the spectral regions 
with our procedure; again, model atmospheres from Kurucz 
(1993) with a grid of different \ogg, T c ff , logiV(Fe) were 
used. 

The results of our "re-determination" of parameters of 
the solar atmosphere are given in Table Q The best fit to 
one spectral region is shown in Fig. From our analysis 
of the solar spectrum we obtained T c ff = 5625 ± 125 K, 
logTV(Fe) = -4.48 ±0.15 dex, V t = 1.2 ± 0.4 km s~\ Here 
and below we used the standard deviation for error esti- 
mates. All these parameters are in good agreement with the 



known parameters of the Sun (Allen 1973) as well as V s = 
Vmacro = 44.3 ± 1.3 km s -1 with initial parameter of con- 
volution. 

4.2 V838 Mon 

The obtained best fits of the synthetic spectra to the ob- 
served spectrum of V838 Mon for February 25 are shown 
in Figs. and Spectral regions containing emission lines 
were removed before we conducted our analysis. Each plot 
shows the fits for a) "the best" T e fl , and b) a case in which 
the effective temperature has been reduced by 1000 K. The 
latter value corresponds more with the T er r obtained from 
photometric measurements (c.f. Munari et al. 2002a). Our 
procedure allows us to choose the best fit from the grid of 
theoretical spectra computed for model atmospheres with 
different log iV(Fe), K , T c g . As we see in Figs.|S|and[7|some 
differences between observed spectra and synthetic best fits 
still remain. In this work we taken into account only the 
iron abundance, abundances of the other elements suggested 
to be solar. Furthermore, lines of ions can not be described 
properly due to too high gravity of used model atmospheres. 

Our main goal was to determine T e ff and logA r (Fe) for 
the dates of our observations. Results of the determination 
of the parameters of the pseudophotosphere of V838 Mon in 
2002 February and March are given in Table [5] 
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- For February 25 we obtained T cff = 5330 ±300 K, 
logiV(Fc) = -4.7 ±0.14 dex and V t = 13. ± 2.8 km s _1 . 

- For the March 26 data the mean values are T e fi = 
4960 ± 270 K, log iV(Fe) = -4.68 ± 0.11 dex, V t = 12.5 ± 
1.7 km s" 1 . 

- And for March 2 the mean values are T cf f = 5540 ± 
190 K, log JV(Fe) = -4.75 ± 0.14 dex, V t = 13.3 ± 3.2 km 

- We obtained V g = 54 ± 3, 47 ± 3 and 42 ± 5 km s" 1 
for February 25, March 2 and March 26, respectively. 

- The / s parameter provides the heliocentric velocity of 
V838 Mon. We obtained V rad iai = -76 ± 3, -70 ± 3 and 
-65 ± 3 km s _1 for February 25, March 2 and March 26, 
respectively. Most probably, we see some reduction in the 
expansion velocity of the envelope. 



5 DISCUSSION 

From a comparison of our results for all three dates we see 
that: 

• The effective temperature for March 26 is somewhat 
lower then for the previous dates. This is an expected re- 
sult, in view of the gradual cooling of envelope. However, 
for March 2 we found a slightly higher value of temperature 
than for February 25. A possible explanation is the heating 
of the pseudophotosphere as result of the third outburst. 

• The microturbulent velocities are very similar and ex- 
tremely high for all three dates. 

• Our analysis shows a lower value of V s for the later 
dates: the effects of expansion and macroturbulence were 
weakened at the later stages of evolution of the pseudopho- 
tosphere of V838 Mon. 

• The iron abundances log iV(Fe) = -4.7 ±0.14 are similar 
for all dates. 

Our estimates of effective temperature are in a good 
agreement with Kipper et al. (2004), although we used dif- 
ferent procedures of analysis. The iron abundance ([Fe/H] = 
—0.4 ) and microturbulent velocity (K = 12 km s _1 ) found 
by Kipper et al. (2004) for March 18 are in agreement with 
our results. 

Our deduced "effective temperatures" as well as those 
in Kipper et al (2004) do not correspond with values ob- 
tained from photometry (T e ft ~ 4200 K). We assume that 
in our analysis we deal with temperatures in the line forming 
region, rather than with the temperatures at photospheric 
levels which determine the spectral energy distribution of 
V838 Mon and the photometric indices. Indeed, the formally 
determined microturbulent velocity Vt = 13 km s~ exceeds 
the sound velocity in the atmosphere (4-5 km s _1 ). This 
means that the region of formation of atomic lines should 
be heated by dissipation of supersonic motions: the temper- 
ature there should be higher than that given in a plane- 
parallel atmosphere of T efi ~ 4200 K. 

Certainly the effect cannot be explained by sphericity 
effects: the temperature gradients in the extended atmo- 
spheres should be steeper (see Mihalas 1978), therefore tem- 
peratures in the line forming regions should be even lower, 
in contradiction with our results. 

Strong deviations from LTE are known to occur during 
the photospheric stages of the evolution of novae and super- 



1.1 



0.9 



0.7 



0.6 



0.5 



0.4 




V 838 Mon 
T eff =5250 K 
T eff =4500 K 



6480 6500 6520 6540 6560 6580 6600 6620 6640 6660 6680 6700 
Wavelength (A) 



1.1 



0.9 



0.8 



0.7 



0.6 




V 838 Mon 



' eff 



=5000 K 
:4500 K 



4 

6300 



6320 6340 6360 6380 6400 6420 6440 6460 6480 6500 
Wavelength (A) 




1.1 

1 

0.9 

X 

! 0.8 
o 
w 

"TO 

E 0.7 
o 
z 

0.6 
0.5 
0.4 

6120 6140 6160 6180 6200 6220 6240 6260 6280 6300 6320 
Wavelength (A) 

Figure 6. The best fits of synthetic spectra to 11-13 orders of 
the observed spectrum of V838 Mon on February 25, found by 
the minimization procedure. 



novae. The main effect there should be caused by deviations 
from LTE in the ionization balance. However, in our case we 
used lines of the neutral iron, which dominate by number. 
We cannot expect a reduction in the density of Fe I atoms 
in the comparatively cool atmosphere of the star. Further- 
more, we exclude from our analysis strong lines with P Cyg 
profiles. Lines of interest in our study have normal profiles. 
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Table 2. Atmospheric parameters for V838 Mon 
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In our work we used model atmospheres computed by 
Kurucz (1993) for the comparative low value of Vt = 2 km 
s _1 . Our estimates of K in the atmosphere of V838 Mon 
are much higher. In general, the temperature structures of 
the model atmospheres of late spectral classes respond to 
changes in K . Nonetheless, in our case, the results depend 
mainly on temperatures in the line-forming region. We can 
use "standard model atmospheres" only in the first approach 
computations. The dependence of our results on Vt should 
be of second order of importance. 

To check this we computed the grid of models with Vt = 
2 and 15 km s' 1 for T off = 4750-5750 K with step 250 K and 
log g = 1. Unfortunately we were not able compute models 
with log g = in case Vt = 15 km s" 1 . A program SAM12 
(Pavlenko 2003) was used. Line opacity treatment was used 
to account absorbtion of atomic and ion lines. Then we carry 
out our fiting procedure for this grid of models for February 
25 data. We found that differences in results less then 0.1 dex 
in iron abudance and 1 km s _1 in Vt . Still this "K factor" 
must be considered in detail in future work based on the 
more sophisticated models. 

In the present work we determined only the iron abun- 
dance. The effects of other element abundances should be 
taken into account as well. However, the obtained effec- 
tive temperatures of V838 Mon are rather low and lines 



of neutral iron are of prime importance. Changes in the 
abundances of other free electrons donors - sodium, potas- 
sium, etc. - cannot affect the ionization equilibrium of iron. 
However, we suggest continued attempts to determine their 
abundances. In addition to their direct importance to the 
determination of the evolutionary status of V838 Mon, they 
are also of importance for the determination of the opacities 
in the stellar atmosphere. Generally, to determine the abun- 
dances of a few elements from modelling of complex blends 
we need to know a) the number and b) relative intensities 
of contributing lines, i.e the abundances of all elements. 



6 CONCLUSIONS 

We have carried out our study in the framework of a com- 
paratively simple model. The static LTE model for the star 
with expanding pseudophotosphere is the rough enough ap- 
proximation. Most of our results are rather qualitative, and 
they should be confirmed and/or refined in the future. How- 
ever, in general we are confident enough in our results, in 
particular: 

• we obtain some evidence of cooling of the line-forming 
region; 
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Figure 7. The best fits of synthetic spectra to orders 14, 16 and 
17 of the observed spectrum of V838 Mon on February 25, found 
by the minimization procedure. 




V 838 Mon 
T eff =5750 K 
T eff =4500 K 



ACKNOWLEDGMENTS 

We thank Drs. Ulisse Munari, Lisa Crause, Tonu Kipper and 
Valentina Klochkova for providing spectra and for discus- 
sions of our results. We thank Dr. Nye Evans for improving 
text of paper. We thank unknown referee for many helpful 
remarks. This work was partially supported by a PPARC 
visitors grants from PPARC and the Royal Society. YPs 
studies are partially supported by a Small Research Grant 
from American Astronomical Society. 

This research has made use of the SIMBAD database, 
operated at CDS, Strasbourg, France. 



REFERENCES 

Allen C.W., 1973, Astrophysical quantities, 3rd edition, The 

Athlone Press, London 
Banerjee D.P.K., Ashok N.M., 2002, A&A, 395, 161 
Bond H.E., et al., 2003, Natur, 422, 405 
Brown N.J., 2002, IAU Circ, 7785, 1 

Crause L.A., Lawson W.A., Kilkenny D., van Wyk F., Marang 

F., Jones A.F., 2003, MNRAS, 341, 785 
Desidera, S., Munari, U., 2002, IAU Circ, 7982, 1 
Evans A., Geballe T.R., Rushton M.T., Smalley B., van Loon 

J.Th., Eyres S.P.S., Tyne V.H., 2003, MNRAS, 343, 1054 
Henden A., Munari U., Schwartz M.B., 2002, IAU Circ, 7859 
Jones H.R.A., Pavlenko Ya., Viti S., Tennyson J., 2002, MNRAS, 

330, 675J 

Kimeswengcr S., Ledcrclc C, Schmcja S., Armsdorfcr B., 2002, 

MNRAS, 336, L43 
Kipper T., et al., 2004, A&A, 416, 1107 

Kolev D., Mikolajcwski M., Tomow T., Iliev I., Osiwala J., 
Nirski J., Galan C, 2002, Collected Papers, Physics (Shu- 
men,Bulgaria: Shumen University Press), 147 

Kupka F., Piskunov N., Ryabchikova T.A., Stempels H.C., Weiss 
W.W., 1999, A&AS, 138, 119 

Kurucz R.L., Furcnlid I., Brault J., Testerman L., 1984, National 
solar obs. - Sunspot, New Mexico 

Kurucz R.L., 1993, CD-ROM 13 

Mihalas D., 1978, Stellar atmospheres, Freeman & Co. 
Munari U., et al., 2002a, A&A, 389, L51 

Munari U., Henden A., Corradi R.M.L, Zwitter T., 2002b, in 
"Classical Nova Explosions", M.Hernanz and J.Jose eds., AIP 
Conf. Scr. 637, 52 

Osiwala J. P., Mikolajewski M., Tomow T., Galan O, Nirski J., 
2003, ASP Conf. Scr., 303, in press 

Pavlenko Y. V., 1997, Astron. Reps, 41, 537 

Pavlenko Ya.V., Jones H.R.A., 2002. A&A, 397, 967 

Pavlenko Ya.V., 2003. Astron. Reps, 47, 59 

Sokcr N, Tylenda R., 2003, ApJ, 582, L105 

Tylenda R., 2004, A&A, 414,223 

Unsold A., 1955 Physik der Sternatmosphcrcn, 2nd cd. Springer. 
Berlin 

Wisnicwski J. P., Morrison N.D., Bjorkman K.S., Miroslmichcnko 
A.S., Gault A.C., Hoffman J.L., Meade M.R., Nctt J.M., 2003, 
ApJ, 588, 486 



• most probably, the line-forming region is heated by su- 
personic motions - our spectroscopic temperatures exceed 
photometrically determined T e ff by ~ 1000 K; 

• we do not find any significant change in the iron abun- 
dance in atmosphere V838 from February 25 to March 26. 

• we derived a moderate deficit of iron logA^(Fe) ~ —4.7 
in the atmosphere of V838 Mon. 
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